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I. INTRODUCTION 



Many of the most frequent hadronic decay channels of B mesons are mediated by the 
quark level transition b — > cdu. The same hadronic dynamics also governs the Cabibbo 
suppressed b — > csu decays. Typical decays of this kind are B — > Dn, B — ► D*n, B — > Dp, 
B -> D*p, B -> Dif, £ -> L>*if, B _-> DK*, B -> D*K*, B -> £> s #~, S -> . . . 

and will be generically referred to as S — > D7T decays. Since these decays are the simplest of 
a complicated array of hadronic channels a great deal of theoretical work has been devoted 
to their understanding 0, Q 3 0, 0, Q, S E H Q 111 EI Q Q • 



After integrating out the W-boson the weak Hamiltonian for B — > Dtx decays is 

G F _ _ 

7i w = ^V cb V* d [C 1 (ii)(d)) V -A{du)v-A + C 2 {p){cib j ) V -A{djUi) v ^ A } , (1) 

where i, j are color indices, and for p b = 5 GeV, Ci(p,b) = 1-072 and C 2 (pb) = —0.169 at 
NLL order in the NDR scheme [3]. For the Cabibbo suppressed TCw we replace d — > s 
and V* d — >■ V* s . It is convenient to categorize the decays into three classes depending 
on the role played by the spectator in the B meson (where "spectator" is a generic term 
for the flavor structure carried by the light degrees of freedom in the B). Class I decays 
receive contributions from graphs where the pion is emitted at the weak vertex (Fig. HIT), 
while in class II decays the spectator quark ends up in the pion (Figs. mCUHl). Finally, 
class III decays receive both types of contributions. Many of these channels have been well 
studied experimentally 17, EM EH HO; 21. 22l|. see Table |TJ Another method to categorize 



these decays makes use of amplitudes corresponding to the different Wick contractions of 
flavor topologies. These can be read off from Fig. ^and are denoted as T (tree), C (color- 
suppressed), and E (^-exchange or weak annihilation). 

Long ago, it was observed that approximating the matrix elements by the factorized 
product {D\{cb)v-A\B) {tt\ (du)v-A\ty gives an accurate prediction for the branching fractions 
of type-I decays, and a fair prediction for type-Ill decays. For all class-I and -II amplitudes 
a similar procedure was proposed |l[. In terms of two phenomenological parameters ai^, 

iA(B° - D+O = aipvr) (D+\(cb) v „ A \B ){7T-\(du) v „ A \0) , (2) 

iA(B° - DV) = ^V cb V: d a 2 (Dn) (n \(db) v „ A \B }(D \(cu) v ^ A \0) . 

Type-Ill amplitudes are related by isospin to linear combinations of type-I and II decays. 
Naive factorization 1 predicts the universal values a\ = G\ + C 2 /N c and a 2 = C 2 + C\/N c . 
Phenomenological analyses testing the validity of the factorization hypothesis have been 
presented in p| , where typically contributions from E are not included. These contributions 
can be modeled using the vacuum insertion approximation which gives the D — ► it form 
factor at a large time-like momentum transfer q 2 = m 2 B . For this reason, they are often 
estimated to be suppressed relative to the T amplitudes by I^ CTl /m 2 0|. 



In this paper we will use the phrase naive factorization to refer to factoring matrix elements of four quark 
operators even though this may not be a justified procedure, and will use the phrase factorization for 
results which follow from a well-defined limit of QCD. 
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FIG. 1: Decay topologies referred to as tree (T), color-suppressed (C), and VF-exchange (E) and 
the corresponding hadronic channels to which they contribute. 



One rigorous method for investigating factorization in these decays is based on the large 
N c limit of QCD. In this limit the amplitudes for type-I decays start at 0(N^ 2 ) while type-II 
decays are suppressed by l/N c (whence the name color-suppressed). The type-I amplitudes 
have a form similar to Eq. (|2J) since non-factorizable diagrams are suppressed, while type-II 
decays simultaneously receive contributions from factorized and non-factorizable diagrams. 
For a typical class-II decay, a Fierz transformation puts the amplitude into the form 

iA(B° - D\°) = ^=V cb V: d { (C a + ^) (D\°\(db)(cu)\B°) (3) 



+2C 1 (D\°\(dT a b)(cT a u)\B )} . 



where the (V — A) ® (V — A) structure is implicit. The two matrix elements have expansions 
in 1/N C which start with terms of order iV c 1//2 and N c 1//2 , respectively 

' -(D^°\(db)(cu)\B°) = F W + + • • • (4) 



Nc Nr 

1 -{D^A{dT%){cT a u)\B") = ±G? + -L^ + 



where ~ iV c , ~ iV c . The Wilson coefficients in Eq. (0) can be assigned scalings 
with N c following from their perturbative expansions C\ ~ 0(1), C*2 ~ which roughly 
corresponds to the hierarchy in their numerical values at The leading terms are the 
matrix elements Fq*\ which factor in terms of large N c form factors and decay constants 

iV 6 1/2 F W ~ (D(* )0 \cu\0)(TT°\db\B°) + (D { * )0 7r°\cu\0)(0\db\B°) , (5) 

plus the matrix elements G^ which are nonfactorizable. The naive factorization assumption 

would keep only F and neglect G±. This approximation is not justified in the 1/N C 

(*) 

expansion since G\ is enhanced by the large Wilson coefficient C\. In either case, no 
prediction is obtained for the ratio of the B — > Dtt and 5 — > D*n amplitudes, 

= A(E° -> D*%°) = (C 2 + Ci/iV c )F * + (2C 1 /N C )G* 1 
A(5Q -> DW) (C 2 + d/iV c )F + {2C 1 /N C )G 1 ' 1 J 
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Heavy quark symmetry does not operate with large N c factorization because for C and E it 
is broken by the allowed exchange of energetic hard gluons between the heavy quarks and 
the quarks in the pion. In contrast, we will show in this paper that expanding about the 
limit E n ^> A this ratio is predicted to be 1 at leading order in A/Q. Here A ~ Aqcd is a 
typical hadronic scale. 

Another rigorous approach to factorization becomes possible in the limit E n ^> Aqcd 
which corresponds to having an energetic light hadron in the final state. In this paper we 
analyze type-II decays using QCD and an expansion in Aqcd/^6, Aqcd/^c, and Aqcd/-Ev 
(or generically Aqqd/Q where Q = {m&, m c , mi, — m c }). We derive a factorization theorem 
and show that E and C appear at the same order in the power counting, and are suppressed 
by Aqcd/Q relative to T. Arguments for the suppression of C by (Aqcd/Q) 1 and E by 
(Aqcd/Q) 1 ' 2 appear in the literature jjfl, but we are unaware of a derivation that is model 
independent. Our leading order result disagrees with the a2-factorization result. Instead the 
amplitudes for B° -> D^°tt° and B° -> p° are determined by the leading light-cone 
wavefunctions (j) WiP , and two new universal B — > distribution functions. Long distance 
contributions also occur at this order in Aq CD /Q, but are shown to be suppressed relative 
to the short distance contributions by an additional a s (Q)/n. 

For type-I decays a color transparency argument given by Bjorken suggested A(B° — > 
D+7T-) ~ (d + C 2 /N c )f 7T F^ D (m 2 7r ) +0(a s (Q)). In Ref. |j it was argued that this factoriza- 
tion is the leading order prediction in the large energy limit E n ^> Aqcd, and in Refs. 0, 0| 
that a s corrections can be rigorously included. This factorization was extended to all orders 
in a s with the proof of a factorization theorem using the soft-collinear effective theory 



■I 

A(B -> £)Wtt) = N w Z(w , fx) I dx T w (x, m c /m b , fx) fa(x, //) + ... , (7) 



o 



where the ellipses denote power suppressed terms. This result is similar to predictions 
obtained from the hard exclusive scattering formalism of Brodsky-Lepage fH], except for 
the presence of the Isgur-Wise function, The normalization factor is given by 2 

AT u) G F V cb V* d / m B \ 

N y ' = — E n f n y/m D (*)m B 1 H • (8) 

The proof of Eq. (J7J) uses the heavy quark limit, so m D = m D * and N = N*. In Eq. (0), 
0tt(x, fx) is the non-perturbative pion light-cone wave function, and £(wo, //) is evaluated at 
maximum recoil v-v' — > wq = (m% + m 2 - )(t) )/(2msm D (.)). The hard coefficient T(*\x,/jl) = 

Cfl R ((4:X — 2)E n , li, rrib), where the ± correspond to the D and D* respectively, and Cfl_ R = 

Cf^ ± C R is the calculable Wilson coefficient of the operators defined in Eq. (}20|) below. 
The renormalization scale dependence of the hard scattering function T(x, fx) cancels the fx 
dependence in the Isgur-Wise function and pion wave function. In this framework there 
is no longer a need to identify by hand a factorization scale. 3 In the language of SCET 



2 Note for longitudinal D*, n-eu* = n-v'. Production of transverse p's is suppressed by A/Q. 

3 In naive factorization the hadronic matrix elements in Eq. (J2J are independent of the scale that separates 
hard and soft physics. The scale dependence in a\ and 02 then causes the physical amplitudes to become 
scale dependent. The parameters a\ and 02 were therefore assumed to be evaluated at a specific scale 
called the "factorization scale" . In other words, the non-factorizablc effects were accounted for by allowing 
ai and C12 to be free parameters that are fit to data. The factorization scale can then be extracted from 
the scale dependence of a\ and 02 0- 
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the scale dependence is understood from the matching and running procedure. 

Eq. (J7J) implies equal rates for B° — > D + n~ and B° — > D* + 7r~ up to the a s (mi,) corrections 
in TW and power corrections. This prediction is in good agreement with the observed data 
for type-I and III decays to tt, p, K and K* as shown in Tables I and II. For two-body type- 
I decays both the large N c and large energy mechanisms make similar phenomenological 
predictions. However, these mechanisms can be distinguished with B — > DX decays where 
X is a multi-hadron state 2 • 



So far, no results of comparable theoretical rigor exist for the color-suppressed type- 
II decays. In fact existing results in B — > Dn and B — > ipK^ do not support naive 



factorization with a universal coefficient 0,2 Furthermore, it has been argued that in 



general factorization will not hold for type-II decays 1711 . 

Using the soft-collinear effective theory (SCET) |25l I2HI]. we prove in this paper a factor- 
ization theorem for color-suppressed (type-II) B — ► DM decays, M = {7r°,p°, . . .}. These 
decays are power suppressed relative to the type I decays, and our results are valid at leading 
nonvanishing order in A/Q. The main results of our paper are 

• The color suppressed (C) and exchange (E) contributions to B° — > D^°ir° are both 
suppressed by A/Q relative to the amplitude (T). The C and E amplitudes are found 
to be of comparable size since the factorization theorem relates them to the same 
perturbative and non-perturbative quantities. Our result is incompatible with the 
naive a-i type factorization. 

• When our result is combined with heavy quark symmetry it predicts the equality of 
the amplitudes for B° -> D°ir° and B° -> D*°ti° (in fact for any DM and D*M). 
This prediction is in good agreement with existing data and will be tested by future 
measurements. 

• Our result gives a new mechanism for generating non-perturbative strong phases for 
exclusive decays within the framework of factorization. For DM and D*M it im- 
plies the equality of the strong phases 5 between isospin amplitudes. Furthermore, 
certain cases with different light mesons M are predicted to also have a universal 
non-perturbative strong phase (f) in their isospin triangle. 

• The power suppressed amplitudes for all color suppressed B — > D^M decays are 
factorizable into two types of terms, which we refer to as short distance (p 2 ~ Em A) 
and long distance (p 2 ~ A 2 ) contributions. The short distance contributions depend 
on complex soft B° — > distribution functions, Sf'^(k + ,£ + ), which depend only 
on the direction of M (the superscripts indicate that two color structures contribute). 
For M = tt,p the long distance contributions vanish at lowest order in a s {Q)/n. 

Combined with Eq. J7J) the results here give a complete leading order description of the 
B — > Dtt isospin triangles. 

In Section |H] we review the current data for B — > Dir decays. The derivation of a 
factorization theorem for the color suppressed channels B° — > D^n and B° — > D^*'°p° is 
carried out in section ITH1 using SCET. Then in section HVI the formalism is applied to decays 
with kaons, B° -> D^°K°, B° -> D^°K*°, B° -> D { * ] ] K~ , and B° -> D { S * ] K*~. In section 
|V]we contrast our results with the large N c limit of QCD and prior theoretical expectations. 
Readers only interested in final results can safely skip sections IlIIl IIVI and|S In section IVT1 
we discuss the phenomenological predictions that follow from our new formalism for color 



5 



suppressed channels. Conclusions are given in IVII1 In Appendix [X] we prove that for tt° 
and p° the long distance contributions are suppressed. Finally in Appendices B and C we 
elaborate on the properties of the jet functions and our new soft B — > distribution 
functions respectively. 



II. DATA 



We start by reviewing existing data on the B — > D^n decays. The branching ratios for 
most of these modes have been measured and the existing results are collected in Table 
Taking into account that the D^ir final state can have isospin / = 1/2,3/2, these decays 
can be parameterized by 2 isospin amplitudes A 1( / 2 , A 3 / 2 : 



A + - 


= A(B° - 


-> D + iT-) 


= ^ A3/2 + \ll Al/2 


= T + E, 


A - 


= A(B~ - 




= V3A 3/2 = T + C, 




A oo 


= A(B° - 




= ^1 A 3/2 - ±A 1/2 





(9) 

Similar expressions can be written for the decay amplitudes of B —>■ D*tt, B — > Dp, B — > D*p 
with well defined helicity of the final state vector mesons. Eq. (JHJ) also gives the alternative 
parameterization of these amplitudes in terms of the amplitudes T, C, E discussed in Sec. HI 
Using the data in Table HI the individual isospin amplitudes Aj and their relative phase 
6 = arg(A 1 / 2 A3 //2 ) can be extracted using 



Br(B -> D^M) = r B T(B -> D W M) = V \A(B -> Z) W M)| 2 . (10) 

Rum- 1 



B pol 

with t bo = 2.343 x 10 12 GeV _1 and r B - = 2.543 x 10 12 GeV~ x . We find 

|Af /2 | = (4.33 ± 0.47) x 10" 7 GeV , 5 Dt = 30.5°+ 7 j , (11) 

\Af /2 \ = (4.45 ±0.17) x 10~ 7 GeV, 

\A?f 2 \ = (4.60 ± 0.36) x 10~ 7 GeV , 5 D * n = 30.2 ± 6.6° , 

\Af f * 2 \ = (4.33 ± 0.19) x 10" 7 GeV. 

The ranges for 5 correspond to la uncertainties for the experimental branching ratios. A 
graphical representation of these results is given in Fig. where we show contour plots 
for the ratios of isospin amplitudes Ri = A 1 / 2 /(v / 2A 3 / 2 ) for both Dix and D*tt final states. 
For B — > Dn an isospin analysis was performed recently by CLEO ljjj including error 
correlations among the decay modes; we used their analysis in quoting errors on 5 Dn . 
For later convenience we define the amplitude ratios 

= A(B° -> D*°M°) jjM/m' = A(B° -> DM°M°) 

A(B° -> D°M°) ' A(B° -> D(*)°M'°) ' 1 ' 

Ai/z _ 3C-E A(B°^D^ + M-)_ C-E 

= Vza 3/2 ~ 1 ~ 2 TTc ' c = a(b- dw°m~) ~ 1 ~ TTc ' 
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Type 


Decay 


XJl^lU J 


|/i| ylU V J 


Decay 


131 l^xU J 


^re v j 


I 


ts 


— > D^ir 


2.68 ± 0.29 a 


5.89 ± 0.32 


B 


— * L> T 7T 


2.76 ± 0.21 


6.05 ± 0.23 


TTT 


B~ 


^D°7T- 




7 7(1 + fl 9Q 


B~ 


-» 7?*°^ 




1 .^iy in u.oo 


II 


B° 




0.292 ± 0.045 6 


1.94 ±0.15 


B° 




0.25 ±0.07 


1.82 ±0.25 


I 


B° 


-> D+p- 


7.8 ±1.4 


10.2 ±0.9 


B° 


-» £>*+p~ 


6.8 ± 1.0 c 


9.08 ±0.68 t 


III 


Br 


-^D°p- 


13.4 ± 1.8 


12.8 ±0.9 


B~ 


-► V 


9.8 ± 1.8 c 


10.5 ±0.97 T 


II 


B° 


^D°p° 


0.29 ±0.11 d 


1.97 ±0.37 


B° 




< 0.56 


< 2.77 



TABLE I: Data on B -> jDWjt and 5 -> DModecays from a Ref. 18], 6 Ref.[ll Q, c Ref.0], 
d Ref. 22], or if not otherwise indicated from Ref.|17l|. ^For B — > D*p the amplitudes for longitudi- 
nally polarized p's are displayed. 



Type 


Decay 


Br(10" 5 ) 


|A| (10~ 7 GeV) 


Decay 


Br(10" 5 ) 


\A\ (W- 7 GeV) 


I 


B° 


-» D+K~ 


20.0 ±6.0 


1.62 ±0.24 


B° 


-> D*+K~ 


20 ± 5 


1.64 ±0.20 


III 


B~ 


— » D°K~ 


37.0 ±6.0 


2.11 ±0.17 


B~ 


-» D*°K- 


36 ± 10 


2.11 ±0.29 


II 


B° 


^D°K° 


5.0+};2 ±0.6 [27] 


0.81 ±0.11 


B° 


-» D*°K° 






I 


B° 


-► D+K*- 


37.0 ± 18.0 


2.24 ±0.54 


B°- 


_> d*+K*~ 


38 ± 15 


2.30 ± 0.45 1 " 


III 


B~ 


-» D°K*- 


61.0 ±23.0 


2.76 ±0.52 


B~ 


-» D*°K*- 


72 ±34 


3.03 ± 0.72t 


II 


B° 


-> D ^* 


4.8±{;J ±0.5 [27] 


0.81 ±0.11 


B° 









TABLE II: Data on Cabibbo suppressed B — * DK^ decays. Unless otherwise indicated, the data 
is taken from Ref.^J. ^ Since no helicity measurements for D*K* are available we show effective 
amplitudes which include contributions from all three helicities. 

where the ratios Ri and R c are denned for each D^M mode. Predictions are obtained for 
the ratios in Eq. ()12j) . including the leading power corrections to Ri and R c . The relation 
Rz — 1 + 0(A/Q) can be represented graphically by a triangle with base normalized to 1 
(see Fig. El in section IVIjl . The two angles adjacent to the base are the strong isospin phase 
8, and another strong phase 0. The usual prediction is that 5 ~ l/Q k fRUlf. and that there 
is no constraint on the strong phase <fi which can be large. In section IV II we show that at 
lowest order the angle <fi is predicted to be the same for all channels in Table EJ and that 5 
can be dominated by a constrained non-perturbative strong phase. From Rj in Eq. (J 12)) we 
note that for a leading order prediction of 5 it is not necessary to know the power corrections 
to the T amplitude. 

A similar analysis can be given for the Cabibbo suppressed B — > DK^ decays. Although 
several of these modes had been seen for some time, it is only recently that some of the 
corresponding class-II decays have been seen by the Belle Collaboration [27] (see Table 
II). For this case the final D^K^ states can have isospins 7 = 0,1, so these decays are 
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parameterized in terms of 2 isospin amplitudes Ar=o,i (f° r given spins of the final particles) 

A 1 = T + C, 

\a x -\a, = C. (13) 

Isospin symmetry implies the amplitude relation among these modes A + ^+Aqq = Aq~, which 
can be used to extract the isospin amplitudes Ao,i and their relative phase S = arg(AoAl). 
Using Gaussian error propagation we obtain 

\A» K \ = (1.45 ±0.62) x 10" 7 GeV, 5 DK = 49.9 ± 9.5° (14) 
\A? K \ = (2.10 ±0.17) x 10~ 7 GeV, 

\A$ K *\ = (1.93 ± 1.49) x 10~ 7 GeV, 5 DK * = 34.9 ± 19.4° (15) 
\A? K *\ = (2.76 ±0.52) x 10~ 7 GeV. 

However, note that scanning the amplitudes A + _, A 00 , A _ in their la allowed regions still 
allows a flat isospin triangle 



A + - 


= A(B° - 


-> D+K~) 


A,- 


= A(B~ - 


-> D° K~) 


A m 


= A(B° - 





III. SOFT-COLLINEAR EFFECTIVE THEORY ANALYSIS 

The key idea of the soft-collinear effective theory |25|, 2(| is to separate perturbative and 
non-perturbative scales directly at the level of operators. The relevant scales have virtualities 
p 2 ~ m 2 ) m 2 c1 E 2 4 (hard), p 2 ~ E M A (intermediate), and p 2 ~ A 2 (soft). The p 2 ~ A 2 scales 
are described by soft (p + ,p~,p L ) ~ (A, A, A) and collinear (p + ,p~,p ± ) ~ (A 2 /E M , E M , A) 
degrees of freedom. We follow the notation in Refs. foL Issl] - 

The weak Hamiltonian Tiw in Eq. (0) is matched onto effective operators containing 
soft and collinear fields. In exclusive processes this matching can be simplified by a two- 
stage procedure [i^. We first match QCD onto a theory SCETi with ultrasoft fields with 
p 2 ~ A 2 , but intermediate-collinear fields with p 2 ~ E^A. This theory gives a simplified 
description of the p 2 ~ E^A exchanges that necessarily mediate interactions between soft 
and collinear particles. Then at a scale fi 2 ~ E^A we match SCETi onto the final theory 
SCETn which has only the propagating long distance soft and collinear particles. This 
procedure determines which factors of a s (/i) belong 4 at the hard scale /z 2 = Q 2 , which 
belong at the intermediate scale fi 2 = EmA, and what non-perturbative matrix elements 
appear. 

Since the collinear fields do not interact with soft fields at lowest order, if one can rear- 
range the fields in the SCETn operator to express it as a product of collinear fields and soft 
fields, the factorization of matrix elements is achieved. This is precisely what happens in 



4 A more accurate statement is that the scale dependence is determined by anomalous dimensions of oper- 
ators in SCETi and SCETn . 



S 



type-I decays, and as we will see also type-II decays, with operators of the form 



Type-I: [h v ,SrS*h v ] [(LW)T' (W^ n )] 

Type-II: [{h v , S)T h v ) {qS)T"{&qj\ [{l n W)V '(W%)] , (16) 



In type-II decays the first and second SCETn operators give short and long distance con- 
tributions respectively We use here the notation in Ref. |26[ so that h v are HQET fields, 
£ n are collinear quark fields, q are soft quark fields, and S, W are soft and collinear Wilson 
lines. Since collinear particles do not connect with the heavy meson states and soft particles 
do not connect with the collinear light meson state, the matrix elements of these operators 
factor into the product of a soft B — > matrix element and a collinear matrix element 
involving M. 

We start by reviewing type-I decays. Using SCET, the factorization of the leading am- 
plitude for type-I decays has been proven in Ref. [9J at leading order in 1/Q (and non- 
perturbatively to all orders in a s ). The operators in Eq. are matched onto effective 
operators at a scale \iq ~ Q 

Y^CiOi^A fdndr 2 [cf (n^Qf (n^+Cf {T X ,T 2 )Qf (n,r 2 )]. (17) 

1,2 j=L,R 

At leading order in SCETj there are four operators [j = L, R] 

Qf (ri,r a ) = [*S?IX>] [(S°W)nr„(Wt6 u) k] , (18) 

The superscript (0, 8) denotes the 1 <8> 1 and T a <S> T a color structures. The Dirac structures 
on the heavy side are r£ R = i/lPl,r with Pr^l — |(1 ± 7s), while on the collinear side we 
have T n = t/iPl/2. The momenta labels are defined by (W'£ n ) W2 = [6(002 — V) W'£ n ]. 
The matching conditions for the Wilson coefficients at tree level at \i — E n are 

Cf fa) = C x + ^ , cf (74) = 2C 2 , 8) (r>) = . (19) 

Matching corrections of order 0(a s ) can be found in Ref. 0]. 

The operators in Eq. (jlKjl are written in terms of collinear fields which do not couple to 
usoft particles at leading order. This was achieved by a decoupling field redefinition |26| 
on the collinear fields £ n — > Y^ n etc. The operators in Eq. (|18|) are then matched onto 
SCETn to give [u i: = r^] 

Qf\u u u; 2 ) = [h^h®] m d) WUT n (W^e ] U] , (20) 

where the collinear and soft Wilson lines W and S are defined in Eq. (jC2|) of Appendix IU1 
At leading order in 1/Q only the operators Q ® R and the leading order collinear and soft La- 
grangians (£ { c \ £i 0) ), contribute to the B~ -> Z>W 7r _ and 5° -> DW+tt - matrix elements. 
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(8) 

The matrix elements of Q L R vanish because they factorize into a product of bilinear matrix 



elements and the octet currents give vanishing contribution between color singlet states 

Note that we take the pion state or interpolating field to be purely collinear and the B and 
D^*> states to be purely soft. Power corrections to these states are included as time ordered 
products. This includes asymmetric configurations containing one soft and one collinear 
quark which involve T-products with subleading Lagrangians |29j. 

Next we consider type-II decays. The matrix elements of the leading order operators 
vanish, (D°ir°\Qf ,8) \B°) = 0. This occurs due to a mismatch between the type of quarks 
produced by Qf' 8 ^ and those required for the light meson state, where we need two collinear 



quarks of the same flavor. The operator Qf' 8 ^ produces collinear quarks with (du) flavor 



Therefore it can not produce a ir° since the leading order SCET Lagrangian only produces or 
annihilates collinear quark pairs of the same flavor. For this reason the leading contributions 
to B° — > D^°7i° are power suppressed. 

In SCETj there are several sources of power suppressed contributions obtained by in- 
cluding higher order four quark operators, higher order contributions from the Lagrangians, 
or both. However, there is only a single type of SCETj operator which contributes to 
B° — > D^°M° decays at leading order. They are given by T-ordered products of the 
leading operators in Eq. (fTSj) with two insertions of the usoft-collinear Lagrangian C^: 

?f 8) = \ Jd A xd% T{Q.f' 8 \o) ,iCf^(x) ,iCf q \y)} . (21) 
Here the subleading Lagrangian is 2^, 3^ 



4? = WUgW ± w)q us - q us (wHgBiW^)(W^ n ) , (22) 

where igp± = [in- D c ,ifl^_]. The two factors of iC^f in Eq. ()21j) are necessary to swap one 
u quark and one d quark from ultrasoft to collinear. In contrast to the tree amplitude, for 
this case both the and operators can contribute. By power counting, the Tj°' 8 ^s 
are suppressed by A 2 = A/Q relative to the leading operators. They will give order A/Q 
contributions in SCETn , in agreement with our earlier statements. 

In Fig. |21 we show graphs contributing to the matching of SCETj operators (a,b) onto 
operators in SCET n (c,d,e). In Figs. |2K D the gluon always has offshellness p 2 ~ E M A due 
to momentum conservation, and is shrunk to a point in SCETn • However, the collinear 
quark propagator in (a,b) can either have p 2 ~ EmA giving rise to the short distance 
SCETn contribution in Fig. |2t, or it can have p 2 ~ A 2 which gives the long distance 
SCETn contribution in Figs. Efc,d. To match onto the short distance contribution in Fig. |2k 
we subtract the SCETn diagrams (c,d): 

(a) + (P) - (c) - (d) = (e) . (23) 

The operators in Figs. |2K,b are from the T-products Tj ' 8 ^ in Eq. (J2TJ), while Figs. I2t,d 

involve the SCET n T-products of in Eq. flZZD, and Fig. |2fe involves Of in Eq. flUJ). 

To generate connected SCETi diagrams from the time-ordered product in Eq. (}2~T]) re- 
quires at least two contractions, of which the minimum basic possibilities can be grouped as 
follows: 
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FIG. 2: Graphs for the tree level matching calculation from SCETt (a,b) onto SCETn (c,d,e). The 
dashed lines are collinear quark propagators and the spring with a line is a collinear gluon. Solid 
lines in (a,b) are ultrasoft and those in (c,d,e) are soft. The (g> denotes an insertion of the weak 
operator, given in Eq. (|18jl for (a,b) and in Eq. (J2UJ) in (c,d). The © in (e) is a 6-quark operator 
from Eq. ()28|) . The two solid dots in (a,b) denote insertions of the mixed usoft-collinear quark 
action £r^ . The boxes denote the SCETn operator cf) in Eq. (|23j) . 

1) Contraction of ffi ffi and the J_ gluon in B^B U ± (C-topology, Fig. Ek), 

2) Contraction of ffl and the i_ gluon in B^B U ± (E-topology, Fig. Eb), 

3) Contraction of ^ ^ and £n £n (topology with two external collinear gluons and 
no external collinear quarks, not shown). 

All more complicated contractions have one of these three as a root. Case 3) only contributes 
for light mesons with an isosinglet component (rj, rf, u, 0), which we will not consider here. 
Each of the SCETt T-products is matched onto SCETn operators at scale \x = /io, and 

[ T L,r] shor t = J iT i dk t duJ k (n, Ho) J m (Ti, kf, UJ k , Ho, fi) 0[°§ (kf, W k , fi) , 

[ T l,r] long = J dk+ duJ i duJ 2 du Cf§ (uji, no) J (0,8) (k + uj, /i , fi) { l^ (u i} k + , u,fi), 

where the subscripts i, £, k run over values 1, 2. Here J, J are jet functions containing effects 
at the p 2 ~ E M A scale and are Wilson coefficients for the SCET n operators O and O. The 
short and Px°r ]iong terms are respectively Fig.Ek and Fig. Et,d (after they are dressed 
with all possible gluons). The fi and fi dependence in Eq. (|2*1|) signifies the scale dependence 
in SCETi and SCETn respectively. The jet functions are generated by the contraction of 
intermediate collinear fields with couplings a s (fio) (where //jj ~ E n A). In general the jet 
functions depend on the large light-cone momenta n coming out of the hard vertex, the 
large light-cone momenta uo^ of the external collinear SCETn fields, and the k^ momenta of 
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FIG. 3: Tree level matching calculation for the L L operators, with (a) the T-product in 
SCETj and (b) the operator in SCETn • Here q, q' are flavor indices and uj ~ A are minus- 
momenta. 



the external soft SCET n fields. No other soft momentum dependence is possible since the 
leading SCETx collinear Lagrangian depends on only n-d us . 

The difference between the time-ordered products T^'^ and the time-ordered products 

Ol'r gi ves the srx quark SCET n operator 0^'^\ whose coefficients are the jet functions 

j(0) 8 ). i n this SCET: -SCET n matching calculation the O^'r graphs subtract long distance 

contributions from the T^'j^ graphs so that j(°' 8 ) are free of infrared singularities. In general 
the matrix elements for color suppressed decays then include both short and long distance 
contributions as displayed in Eq. ()24|). However, for the isotriplet it and p a dramatic 



simplification occurs at leading order in Cj\ In this case it can be proven that the long 

distance contributions [Tj ]i ong vanish to all orders in the a s couplings in SCETi , and 
with the a s couplings in SCETn treated non-perturbatively. The proof of this fact uses 
the G-parity invariance of QCD and is carried out in Appendix [XJ At leading order in 
the coefficients C^r the M = it, p factorization theorem is therefore more predictive since 

possible long-distance contribution from Oj are absent. Most of the following discussion 

will focus on Oj l \ but is fully included in the final factorization theorem. 

In the SCETn diagrams in Fig. |2t,d a power suppressed four quark Lagrangian appears. 
It is similar to an operator introduced in Re f. |3ll| . and can be obtained from T{iC^,iC^} 
in SCETi by a simple matching calculation |32|. Summing over flavors q, q' we find 



j=L,R uj k+ 

Lf{u>Mx) = Y^i&Whtt^tiFU [{Q'S) k+ iP,{^q) k+ ](x). (25) 



g,q 



In Eq. (|2~5|) the soft momenta labels are defined by (S^q)k+ = [5{k + —n-P) S^q], and the po- 
sitions (x + ,x~,xj_) ~ (1/A, Q/A 2 , 1/A). For the soft fields the x~ coordinates encode small 
residual plus-momenta, and for the collinear fields the x + coordinates encode small residual 
minus-momenta. Thus, we used the summation/integration notation for label/residual mo- 
menta from Ref. 33||. The operator Lf\uj, k + , x) has the same form as Eq. (|25|) except with 



color structure T a ®T a . At tree level the coefficient functions are given by the calculation 
in Fig. El 

7 (o) M+) = _Cf_ Ana s (p) (8) = J_ Ana s (p) 

v ; 2A C uok+ ' v ; 2A C uok+ V ; 
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Beyond tree level they obtain contributions from loop diagrams with additional vertices. 

In terms of the operator in Eq. the SCETn operators that contribute to [T- ]i ong in the 
factorization theorem are 

0f' 8 \uj h fc+ v, n) = Jd 4 xT Qf' 8 \uj u x = 0) iL^\u, fc+ x) . (27) 

The operators O generate the diagrams (c) and (d) in Fig. |21 

At any order in perturbation theory the jet functions J from the C-topology and E- 
topology generate one spin structure, and two color structures for the SCETn operators. 
For the six quark operators we find 



Of\ki^ k )= h^Y)h^{dS) kt iP L {S^u) kt (LWUT c (W^nU , (28) 



;(<o 



of\kt,u k ) = [(h^s)r^T a (s^h^) {dS) k +iP L T a {tfu) kt \ [(Lwur^^u 

where here the d, u, hff , and hv fields are soft, and the fields are collinear isospin 
doublets, (&\& d) ). In Eq. ®®T h LR = iP L . Ft as in Eq. (|18|) . while for the collinear isospin 
triplet T c = r 3 ^P^/2. 5 We do not list operators with a T a next to T c since they will give 
vanishing contribution in the collinear matrix element. For light vector mesons the spin 
structure T c only produces the longitudinal polarization. This result follows from the quark 
helicity symmetry of JC^ and is discussed in further detail in Appendix [B] 

In position space the O^ 1 are bilocal operators, with the two soft light quarks aligned on 
the n M light cone direction (x _ = \n^n ■ x, y~ = ^n^n ■ y) passing through the point x = 

(h { :)S)T h (S^h^) (dS) r+ T q (S^u) e+ = (29) 

The gluon interactions contained in matrix elements of Oj 0,8 ' ) include attachments to the 
light quarks q, to the heavy quarks h v y, and to the Wilson lines S n as shown in Fig.HJ The 
interactions with h v y have been drawn as Wilson lines S v y along v , v ' j34|]. Even though 
we have factored the collinear and soft degrees of freedom in the two final state hadrons, the 
presence of the soft Wilson lines bring in information about the vector This allows the 
soft operators oj to be non-trivial functions of n • kj, n • v, and n ■ v', and this information 

gives rise to a complex phase in the soft functions S^'^ as shown in Appendix O Thus, 
the S n Wilson lines are directly responsible for producing final state interactions, and the 
soft fields in oj 0,8 ' ) encode non-perturbative rescattering information. 6 This makes good 
sense given that the soft gluons in the 5„'s were originally generated by integrating out 
attachments to the collinear quarks and gluons making up the light energetic hadron. 

The above procedure provides a new mechanism for generating non-perturbative strong 
phases for exclusive decays within factorization. In the soft B — > matrix elements 



5 There are also isosinglet contributions with r c = ijiP^JI. 

6 Note that in semi-inclusive processes a different mechanism is responsible for the phases in single-spin 
asymmetries which has to do with the boundary conditions on Wilson lines j3^|. 
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FIG. 4: Non-perturbative structure of the soft operators in Eq. (|29|) which arise from of' 8) . Wilson 
lines are shown for the paths S n (x, 0), S n (0, y), S v (—oo, 0) and S v >(0, oo), plus two interacting QCD 
quark fields inserted at the locations x and y. The S v and S v > Wilson lines are from interactions 
with the fields h v and h v > fields, respectively. The non-perturbative structure of soft fields in O^f^ 
is similar except that we separate the single and double Wilson lines by an amount x±. 



the information about the light energetic meson is limited to its direction of motion n/\ 
Since these matrix elements know nothing further about the nature of the light meson, 
these strong phases are universal. In particular the same strong phase <fr is generated for the 
decays B — > D^ix and B — > D^p. (We caution that this is not the isospin strong phase, but 
rather a different angle in the triangle.) The same mechanism produces another universal 
strong phase for color suppressed decays to DR(*>°, and a third for decays to D S K^~. The 
different phases in the three classes arise in part due to the appearance of different moments 
of the matrix elements of the soft operators. However, for the kaons there are additional 
long distance contributions to the strong phases from [T]i ong , which make the universality of 
the phase from [T] short hard to test. A more complete set of phenomenological predictions 
is given in Section IVIl including a comparison with existing data. Further details on the 
properties of the soft functions S^ 0,8 ) are given in the Appendix IU1 

The matrix elements of the short distance operators in Eq. (|28p factor into products 
of soft and collinear parts, respectively. The collinear part of the matrix elements are simply 
given in terms of the light cone wave function of the light meson. For the it and p the 
definitions are [we suppress pre-factors of f Q dx5(u>i — xn-p M ) <5(u; 2 + (1 — x)n-p M ) on the 
RHS] 7 

KlfaWUfaTsiW^JM = -iV2Ufi-p v <p v (ji,x) , (30) 
(pl{£)\{lnWU$Tz{W^ n )^\V>) = iV2f p m p n-e*cf> p (v,x) 

= i\/2 f p n-p p 4> p (p,x) . 

In the last equality we have used the fact that at this order the collinear operator only 
produces longitudinal p's, for which m p n-e* L = n-p p . 

Since it no longer contains couplings to energetic gluons, the soft part of the matrix 
elements of oj 0,8 -* can be constrained using heavy quark symmetry. In other words, heavy 



7 Our vector meson states are defined with an extra minus sign relative to the standard convention. 
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quark symmetry relations can be derived for matrix elements of soft fields. The constraints 
can be implemented most compactly using the trace formalism of the HQET (36[. First 
consider the matrix element of the soft fields in Oj 0,8 ^. For we have 

i(^)r (s^) {dS) kt iP L (&u) k} \B°(v)) _ (c) 



^/m B m D 



Tr [H^/TH^X^] , (31) 



where X^ ' = X^(kj~, n, v , v ') and we use the standard relativistic normalization for the 
states (and note that the LHS is independent of m^ c in the heavy quark limit). An identical 
equation holds for Oj with an X^ 8 \ In writing the trace formula in Eq. (|31|) we have used 
the fact that the d and u quarks must end up in the B and states. 8 The heavy mesons 
(D,D*) and (B,B*) are grouped together into superfields (3a |. defined as 



H v = ^l(P:^ + P vl5 ). (32) 

Now X(°' 8 ) are the most general structures compatible with the symmetries of QCD. They 
involve 4 functions af'^lkf , k^, v -v', n-v,n-v') 

X^ = a^iP L + a^iP R + af 8) P L + af> 8) P R , (33) 

Structures proportional to j> and f can be eliminated by using H v f = —H v , etc. 

The presence of four functions in Eq. (|3^j) would appear to restrict the predictive power 
of heavy quark symmetry. However, using the properties of H v and H v i and the fact that 
the two-body kinematics relates n to v and v' via tubV = mov' + Emu, it is easy to 
see that the four functions dj appear only in two distinct combinations. (Note that we 
are taking m M /m B ~ A/m B <C 1.) For T*l R they give soft functions Sl,r defined as 
Sl = (n ■ f')( a i — a 3/2) — 0,4/2, Sr= [n ■ v')(a 2 — a 4 /2) — a 3 /2 and 



(^(^[(^^^(^/.(^(^^^(^^l^^)) i(0) 



(D- o (^^ g )l(^ ) g)^(gt^))(Jg) fc ^p L (^ )fc2+ | J BO (^; ) ) , 

where the ± for the D* refers to the choice of Pl or P B . Identical definitions hold for the 
matrix elements of the color-octet operators which give S^ R (k^). We will see in section IVT1 
that the result in Eq. ()34|) relates decay amplitudes and strong phases for B° — > D°M° 
and B° — > D*°M° at leading order in the power expansion, and up to terms suppressed by 
a s {Q)/-K. If one takes n ■ v — 1, then n-v' = ms/m^ , v-v' — (m B + m 2 DM ) / (2mBm D (*)) . 
The D,D* variables are equal in the heavy quark limit. 

For the long distance operators Oj the same set of arguments in Eqs. (jSIHSU) can be 
applied except that now we must add terms cl^'^^Pl + cl^' 8 ^$±Pr to X^ 0,8 \ and the dj's 



8 The matrix element of the analogous soft operators with (uu)+(dd) would contain a second term in Eq. (jHlJI 
of the form Tr [77;, TH^ X}Tt [Y], which arises from contracting the light quarks in the operator. These 
types of traces also show up for power corrections to B° — > D^*' + M~ and B~ — > D'*'°M" . 
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can also depend on x\. The functions analogous to S^' R are defined as jj (A; + , x±, s* D *). 
In this case the D and D* decompositions are no longer related since the matrix element 
involves both n-e* and x±_ -e* terms for the D*. Thus, due to the long distance contributions 
for light vector meson we must restrict ourselves to the longitudinal polarization in order 
to have equality for the D and D* amplitudes. In the case of the p this restriction is not 
important since the long distance contributions vanish (see Appendix However this 
observation does have phenomenological implications for decays to K*'s. 

We are now in a position to write down the most general factorized result for the amplitude 
for the decays B° — > D^°M°. Combining all the factors, this formula contains the soft 
functions S^ 0,8 '(kf, kf) from Eq. (|34jl . the jet functions from Eq. 1)240. and the Wilson 
coefficients Cf§ from Eq. (JTZJ). In J®(n,kf ,Uk) we can pull out a factor of 8(ti — t 2 — 
uji + 002) by momentum conservation. This leaves the variables n + t 2 = 2Em(2z — 1) and 
UJ1+U2 = 2Em{2x~1) unconstrained, which give convolutions with the momentum fractions 
z and x respectively. In defining J^(z, x, kf) we multiply JW(rj, kf, Uk) by cji — 0J2 = n-pM- 
All together the result for the B° — > Z)(* )0 M amplitude is 

A£ W = iV A/ J dx dzjdkfdkf Cf{z) J {i \z, x, kf, kf) S<£\kf, kf) 4> M (x) (35) 

±cff(z) J^(z,x,k+,kf) S$(k+,kf) <f> M (x) 

"r^-long ) 

where we sum over i — 0, 8 and the po, P dependence is as in Eq. (j21j). The A^ ( ^ M in 
Eq. 1)3*5)) denotes the contributions from the matrix elements of the SCETn time-ordered 
products [T]i 011 g. Also the ± refer to D/D*, C^ R (z) = C^' r (ti + T2, Em, rrib, mc> /•*)> an d 

n m = G F v«y: d h[ ^ mBmDM (36) 

The normalization factor is common since mo = m* D and n-e^ D *^ = n-v ' . This follows since 
the M's produced by Oj 0,8 ^ are longitudinally polarized. 

The long distance amplitudes also obey a factorization theorem which can be derived 

by examining the matrix elements of the 0^' R operators in Eq. (|27|). First factorize the 
collinear fields into the matrix element with the M and the soft fields into the matrix element 
with the B, £)". The independence of the collinear propagators on the residual soft minus- 
momenta leads to a S(x + ) and the independence of the soft propagators on the residual 
collinear plus- momenta leads to a S(x~) (somewhat similar to the calculation for B — > X s 7 
as described in Ref. [26]). The result is 

Aong M = K Jdzjdk+dLojd^Affiz) J«(^ + ) ^f{k + ,x ± ,e* D ,)^{z,u;,x ± ,e M ) 

±C®(z) J (i) (^ + ) <S> { £(k + ,x ± ,e* D *)^(z,u;,x ± ,e* M )] . (37) 

where the ± is for D and D* and we defined the non-perturbative functions in a way which 
gives the same prefactor as in Eq. 023). HereCf^ are the Wilson coefficients of the weak op- 
erators in Eq. ()2U|). and the jet functions t /°' 8 ' ) are the coefficients of the SCETn Lagrangian 
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in Eq. ([25)1 . The B and are soft and collinear matrix elements from the operators O 
and are given by [with prefactor dz8{uj\ — zti-pm)5(uj 2 + (1 — z)n-pM) for \1> 



(0)! 
M . 



M°(pm, e M ) [(^WO* JPl(^(«))^J (0 ± ) [(^ ) ^)^Pl(^W )„] (x ± ) 
= if M /V2^)(z,uj,x ± ,e* M ) , 

DMV,^o|[(^)r^^ b ( 

= ±y/m B m DW $^ R (k + , x ± , e* D .) , 



(38) 



and ^f\J and ij are defined by analagous equations with color structure T a ® T a . The ± 

is for Pl and Pr respectively. In a more traditional language the A®^ M contributions might 
be referred to as "non-factorizable" since they involve an direct xj_ convolution between non- 
perturbative functions. Eqs. ()35|) and (J3~7j) are the main results of our paper. Additional 
details about the derivation of Eq. (|3~7jl will be presented in Ref. |37| |. 

Using the SCETn power counting in rj = A/Q we can verify that the short and long dis- 
tance contributions to the factorization theorem are indeed the same order. The coefficients 



77°. The results in Eqs. (|26j) and (|42|) for the jet functions imply JW ~ 1/A 2 and 

J y "' ~ 1/(QA). Furthermore, 0m ~ ^° from the definitions in Eq. fl3()l) . For the soft function 
in Eq. (|34j) we get (77 s / 2 ) 4 from the fields, r/ -3 from the states, times r]~ 2 from the delta func- 
tions indicated by the momentum subscripts. This gives S(kf, k£) ~ rj, ie. S(kf, k£) ~ A. 
A similar calculation for the collinear and soft long distance matrix elements in Eq. (|38|) 
gives ^m' 8 '' ~ A 2 /Q an d ^i'R ~ A. In the factorization theorem the measures have scaling 
(dk^dk^) ~ A 2 and (dk + d 2 x±) ~ 1/A. Combining all the factors for the short distance 
amplitude gives (A)(A 2 )(1/A 2 )(A)(A°) = A 2 , while for the long distance amplitude we find 
(A) (1/A) (1/A) (A) (A 2 ) = A 2 also. Therefore, both terms in A(£ 

are the same order in the 
power counting. They also give the complete set of contributions at this order. 



For numerical results with M = it, p the A® * M contributions are very small since taking 



C^' R (z) independent of z gives A^' M = as shown in Appendix |XJ This implies that 



A B Q * g M / A m ~ a s (Q)/ir, and together with the helicity structure of the jet function dis- 
cussed in Appendix |B] implies that the production of transverse p mesons is suppressed. In 
Section [VTI we explore further phenomenological implications. 

Next tree level results are presented for the jet functions J(°' 8 \ The SCETi graphs in 



Fig.Elare computed with insertions of Qj ' and taking momenta —k\ and — ki for the initial 
and final light soft antiquarks, together with momenta p% and P2 for the collinear quark and 
antiquark. The diagrams in Fig. Ek,b with insertions of {Qj , Qj } are 



C: 9< 



P L {1, T B }u^) (uj d) lu P L j/2{l, T B }T A ^ ± vj u) ) (^TY^ 1 ) 
[n- (ki—k 2 )+ie] [n-p 2 n-ki+ie] 



[n- (ki — k 2 )+ie)[—n-px n-k 2 + ie] 



Adding these contributions with factors of and Cf 1 to distinguish the two color struc- 
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tures, and then Fierzing gives 



(o) 



u$iP L u® vi d) iP L vi u A (40) 



X- 1 



N c \ [n- {k\— k 2 ) + ie] \n-p 2 n-k\ + ie] [n- (k± — k 2 ) + ie] [— n-p\n- k 2 + ie] 



(8) 



v$iP L T a u® vfHP L T a v^ 



I -(d) -ir) (A -i u )-iT-> ( u ) 

x TTa s ( U n TjiPlVn Un fPhVn 



iV 2 y [n- (k\ — k 2 ) + ie][n-p 2 n-k\ + ie] [n-(ki — k 2 ) + ie][— n-p\n-k 2 + ie] 

where Cp = (iV 2 — 1)/(2N C ) and we set C^' 8 ^ = 0. The first term in each round bracket 
originates from the C— type graph (Fig. |2K) and the second term from the E— type graph 
(Fig. |2t>). It is convenient to group the result into isosinglet and isotriplet terms for the 
collinear spinors. Since the ir° and p° have definite charge conjugation we can freely inter- 
change the positive momenta n • p\ «-> n ■ p 2 , so a factor of 1/n • p± can be pulled out front. 
For the terms in round brackets we find 



1 [u^^P L v n d) - u^P L vi u) ] 1 [[u^PA d) + u^^P L v^} n-(k 2 + k, 



[n-ki + ie][— n-k 2 + ie] 2 [n- (ki~ k 2 ) + ie] [n- k\ + ie] [— n- k 2 + ie] 



(41) 



For B° — > D^n and B° — > D^*'°p° where we have isotriplet M°'s the contributions 
from the SCETn diagrams in Figs. Efc,d cancel. Thus, the denominator in Eq. (|41|) directly 
gives the tree level isotriplet jet functions 

•/<»>(;, *f , k } ) = J^M^L §^ T — ] , (42) 

N c x[n-ki + ie][—n-k 2 + ie] 



t( 8 )( 7 t h + h+) 



27ra s (yu) <5(z — x) 



A?" 2 x [n- k\ + ie][— n-&2 + ie] 



where n-p\ = x n-pM- These jet functions are non-singular given that the non-perturbative 
soft function S(kf ,k£) vanishes for k^ = or k 2 = 0, and that 4> n ,p(x) vanishes at x = 
and x = 1. On the other hand for isosinglet M°'s the result in Eq. ()41|) has a singular 
denominator 1 / [n-(ki — k 2 ) +ie] . The singularity occurs when the collinear quark propagators 
in Figs. |2K,b get too close to their mass shells, ie. when n-(ki — k 2 ) < A 2 / Q. This singularity 
is exactly what is canceled by subtracting the SCETn diagrams in Figs. I2t,d, which then 
gives a non-singular isosinglet jet function. 

Next we consider the result for the factorization theorem for M = it, p with these tree 
level jet functions. Taking the matrix elements of the O^' 8 ^ operators, the collinear part 
factors from the soft operators as explained above. Their matrix elements are given in terms 
of the M° light cone wave function, and the S^°' 8 \k + , l + ) functions, respectively. This gives 
the explicit result for the B° — ► D^ ir° and B° — ► D^°p° decay amplitudes, at lowest order 
in the matching for C and J 



A(B° - 


+ £>W°7r°) 


= N£ 


A(B° - 


- £><*V) 


= < 



^(y s (p )C F ( o) (0) 27ra s (p ) (8) (8) 1 x 



iV c ' AT? 



(43) 
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We choose to evaluate Cf' 8 \ s^ 0,8 \ and (x^ 1 ) at the common scales /x = /i ~ VEJV since 
one of the hard scales m 2 c is not much different than E n A. In Eq. (J4"3|) the convolutions of 
the soft and collinear matrix elements are defined by 

s (o l8) = | s (o )8) | e #(o.B) = r dk + dk + s<£*\kt,i4,n) 



(k+ + ie)(-k+ + le) ' 
(•'" — 1 )m = A'- • (44) 



o 



^From Eq. (|44j) we can immediately verify the result of the power counting for operators 
described earlier. Since {x~ 1 )m ~ {x°) m ~ A , comparing Eqs. (|7|8jl and (|4Hjl we see that 

— t=- v ~ 47ra s (/i ) — t4 — ~ 47ra s (/in) - — ~ 47ra s (/in) QCD , (45) 

where we have used the standard HQET power counting for the soft matrix elements to 
determine that s^ ' 8 ) ~ A QCD . Thus, the ratio of type-II to type-I amplitudes scales as 
A/Q just as predicted. Due to the factor of 4tt the suppression by a s does not have much 
effect numerically. The 4ir arises because the a s is generated at tree level. It is expected 
that perturbative corrections to the matching for C and J will be suppressed by factors of 
a s (Q)/7r and a s ^E n A)/n respectively. In Eq. (|4"H|) grouping g 2 N c ~ 1 gives an extra factor 
of 1/N C , so with this counting the ratio is color suppressed as expected. 



IV. ADDING STRANGE QUARKS 



In this section we consider how the factorization theorem derived in section ITTT1 is modified 
in the case of color suppressed decays involving kaons, which include B° — > D*K~, B° — > 
Di* ] K*-, as well as the Cabbibo suppressed decays B° -> D^°K° and B° -> D^°K*°. 

If strange quarks are included in the final state then operators with different flavor struc- 
ture appear. In the exchange topology we can have the production of an ss pair (as shown 
by the s-quarks in brackets in Fig. 12b). This gives SCETn six-quark operators 



of {#,«>„) = (/^rjT^s^f ) {dS) k+i iP L T\sh) kt \ ki s) wur c (w^e 



(46) 



which mediate B° — ► D . For the long distance contribution we take flavors q' = d 

and q = s in the Lagrangian in Eq. ([25)1. which leads to s,s quarks replacing u, u quarks 

in % LR . The result for the factorization theorem is then identical to Eqs. (J35J) and (pTTj). 
except that only the E-topology contributes. For this case the long distance contribution is 
not suppressed, and serves to regulate the singularity when matching onto the E-topology 
jet functions j(°' 8 ) = J^ 0,8 -*. Further discussion of the singularities is left to Ref. j3?|]. The 
hard coefficients Cr°'J are the same as in the previous section. 



L,R 



The remaining difference for B° - 
light-cone wavefunctions for K~ , K°, 



are the non-perturbative functions. The 

K*~ , and K*° are [with q = u,d, u>% — n-px s , u 2 = 
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—n-px q , and a pref actor as in Eq. (JHUj)] 



-2if K n-p K <j> K (ji,,x a ) , 
-2if K * m K * n-e* 4> K *(fi,x s 
-2if K * n-p K * (j) K *{n,x s ) . 



(47) 



The collinear functions \E f ^' 8 ' ) also depend on the light meson M. The non-perturbative 



soft functions involve strange quarks and are also different from section ITTT1 5 



U0,8) 



§l R and 

®l'r ~^ ®l'r ■ The non-perturbative functions are related to those in the previous section 
in the 577(3) flavor symmetry limit. However, the jet functions are not related in this limit, 
they differ since different topologies contribute. This leads to different convolutions over the 
non-perturbative functions. 

Next consider the Cabibbo suppressed b — > csu transition with the color suppressed 
topology (as shown by the brackets in Fig. |2k). For the six quark operators we have 9 



1^2 



(48) 



of\kt,u k ) = (h$s)r*T a (sihV)(dS) kt jtP L T%siu) kt (^wur c (w^. 



;(<=) 



which mediate the decays B° — > In this case the SCETn Lagrangian in Eq. (J25|) 

has the same flavor structure as in section ITTT1 Since only the C-topology contributes the long 
distance contribution is not suppressed in the factorization theorem, and the jet function 
j(o,8) J^ ,8 \ For both the short and long distance non-perturbative functions the change 
of flavor appears only through the collinear quarks in the weak operator, so the collinear 



functions depend on the but the soft functions S^'r and are identical to those in 



,(0,8) 



section ITTT1 (However, now J^' 8 ^ appears, so the moments over the soft function 5^ ; will 
be different.) Finally note that if we allow a strange quark in the initial state (for T^-decays) 
then the E-topology can also contribute and more operators are generated. 

Due to the non-negligible long distance contributions the number of model independent 
phenomenological predictions for kaons are more limited. The main predictions are the 
equality of branching fractions and strong phase shifts for decays to D versus D*. For 
M = K°, K~ an identical proof to the one for 7r° and p° can be used. For the vector mesons 
the proof can also be used if we restrict our attention to longitudinal polarizations, so the 
final states D^Kf are related, and so are D^*'KT, . The factorization theorem allows for 
transversely polarized kaons at the same order in the power counting, but only through the 
long distance contribution. 



L,R 



^(0,8) 



V. DISCUSSION AND COMPARISON WITH THE LARGE N c LIMIT 

It is instructive to compare the N c scaling of the different terms in the SCET result 
Eq. (|35|) (or Eq. (J43|) ) with that expected from QCD before expanding in 1/Q given in 



9 Note that the flavor structure was not distinguished in naming the operators in Eqs. Ij28l46l48|l . This 
should not cause confusion since they always contribute to different decays. 
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Eq. ©. Combining the matrix elements in Eq. © written in a form similar to Eq. ()43|) 
gives the decay amplitude at leading order in 1/Q as 



A(B C 



D°M°" 



N, 



M 



Ci + 



9i 



N, 



(F + 2G 1 ) + - 



(49) 



+ iV M C 2 



F + j^(-F + F 2 -2G 1 ) + 



The ellipses denote power suppressed terms. This reproduces the 1/N C expansion of the 
SCET amplitude in Eq. (J43)) with the identification 







G\ = —7ra s CF s 



(0)| 



F 2 - 2G?i = 2na s s (8) 



(50) 



where s(°> ~ N® and ~ AT C . This implies that the factorizable term F is power sup- 
pressed in the limit of an energetic pion relative to the leading order amplitude in Eq. (|43|). 

The naive factorization approach in Eq. (J2J) keeps only the F term, which is expressed in 
terms of the B — > 7r form factor in the large iV c limit. We comment here on the form of this 
contribution in the effective theory. They appear in the matching of the (db)v-A(cu)v-A 
operator onto SCETi T— products such as 



(i) 



T. 



(4) 



(51) 



where the operators Q( 2 ' 3 ) contain one usoft light quark. From the leading order opera- 
tors in Eq. (J2Uj) they can be constructed by switching £ n — > q to give Q^ 2 \ and adding a 
further WHD ± W to get Q (3) . Their precise form is different depending on whether they 
are introduced by matching from the color-suppressed (C) or the W -exchange (E) graph. 
Schematically 



C-type: 



E-type: 



Off 

qS 2) 



1 



(52) 



[(^|^)._. pf 



hU\ 



— [^WUY c h^h^[W^p^ c WUY h u] 

Tflf. 



Q 



(3) 
la 



o (3) 

^16 



[rfr^][^ ) r c ^(wt^ ±c | e w )a;] 

— [dT h h^] [WUp ±c WUT c (W^^ 



(53) 



The presence of the usoft quark field g in these operators introduces an additional sup- 
pression factor of A 2 , such that the T— products T 12 are 0(A 4 ) ~ A 2 /Q 2 down relative to the 
operators Q[ 5 in Eq. (Note that since the form factors enter as time ordered products 



we do not expect a different a s suppression for lf$ relative to those in Eq. (EH) 0.) This 
explains the absence of the Fq contributions at order A/Q, as noted in ([50)1 . Although Fq is 
part of the leading order result in the large N c limit, it is sub leading in the 1/Q expansion. 
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After soft-collinear factorization, the T-products (|51|) match onto factorizable operators 
in SCETu . For example, the C— type time-ordered product containing Q\J gives (schemat- 
ically) 

T 2 (4) - f du x du 2 J{u h fc+) \(dS) kl T(S^h^)} [h^r h u]{(^WUT c (W^^U] (54) 



Apart from the (cu) soft bilinear, this is similar to a factorizable operator contributing to 
the B — > 7r form factor [29]. The presence of the D meson in the final state implies that the 
matrix element of the soft operator in Eq. (J54")l is different from that appearing in B — > n. 
Therefore, naive factorization of type-II decay amplitudes, as written in Eq. (J2J), does not 
follow in general from the large energy limit. Still, in the large N c limit, the matrix element 
of above can be indeed expressed in terms of the B — > ir form factor, as required by 
Eq. © 

Recently an analysis of color-suppressed decays was performed using the "pQCD" ap- 
proach working at leading order in an expansion in m D (*)/m B and Aqcd/^dw [38[. This 
differs from the expansion used here in that we do not expand in m D {*) /mg. The non- 
perturbative functions in their proposed factorization formula include the light-cone wave- 
functions <p^\x3), 4>d{x2) and a B light-cone wavefunction that depends on a transverse 
coordinate </>b(xi, 61). This differs from our result which involves a B — > D function 
Sik^k^) and also has additional long distance contributions, A®* g M , at the same or- 
der in our power counting. Our long distance contributions are "non-factorizable" in the 
sense that the non-perturbative functions $^ R (fc + , x±) and ^/^(z,u,x±) communicate di- 
rectly through their x± dependence without going through a hard kernel. In Ref. [38^ 
strong phases only occur from the perturbative /1q ~ Em A scale, whereas we also find non- 
perturbative strong phases from the A 2 scale (in S(kf , k^)). The non-perturbative phases 
are expected to dominate in our result. Finally, the results in Ref. |38] do not manifestly 
predict the equality of the D and D* amplitudes since at the order they are working contri- 
butions from different B — ► M form factors show up. For example their pQCD prediction 
Br(B° -> D*°p°)/Br(B° -> D°p°) = 2.7 is much different than the prediction of 1.0 that 
we obtain in the next section using heavy quark symmetry. 

The time ordered products presented in this paper in Eq. (J21|) are only A/Q down from the 
class-I T amplitudes. Therefore, they give the dominant contribution to the color-suppressed 
and W-exchange amplitudes in the limit of an energetic pion (A/Q 1). This is a new 
result, not noticed previously in the literature. The power counting of "factorizable" Fq 
type contributions are indeed suppressed by A 2 /Q 2 in our analysis in agreement with the 
literature. However, these terms do not give the dominant contribution. 



VI. PHENOMENOLOGICAL PREDICTIONS 

A factorization theorem for color-suppressed B° — > D°M° decays was proven in Sec- 
tion Unl and extended to decays to kaons in Section ITVl The amplitudes at leading order in 



22 



Aqcd/Q with Q = {mb,m c , E n } have the form 
A 00 = A(B° -> D«°M ) 

= N M j'dxdzjdk+dkt [cf{z)Sf{k+)±C%{z)S [ £{k+)\ J®(z,x,k+)4>M(x) 



i=0,l 

+^Sng M > (55) 

where the sign ± corresponds to a D° or D*° meson in the final state, respectively. In this 
section the implications of Eq. (J55j) for the phenomenology of color suppressed decays are 
discussed. One class of predictions follow without any assumptions about the form of J: 

• Heavy quark symmetry relates the nonperturbative soft matrix elements appearing in 
the B° — > D°M° and B° — ► D*°M° decays with the same light meson at leading order 
in a s (Q)/ir. This implies relations among their branching fractions and equal strong 
phases in their isospin triangles. 

These relations are encoded in the ratios Rq 1 in Eq. (|12p. A second class of predictions 
depend on using a perturbative expansion of J in a s (fio) for /1q ~ Em A: 

• Using a perturbative description of J the amplitudes and strong phases for decays to 
different light mesons M can be related at leading order in a s (fj, )/ir. 

These predictions are encoded in the ratios R^ M , R c , and strong phase in Rj, as defined 
in Eq. (|12|). We consider the two classes of predictions in turn. 

First, consider relations between color-suppressed B — > DM and B — > D*M decays with 
the same light meson. At tree level in the matching at the hard scale fi ~ Q, two of the 
Wilson coefficients vanish C^' 8 ^ = 0. Therefore both amplitudes for D and D* contain only 
the soft functions S^' 8 ^ (A^*) appearing in the same linear combination. This implies model- 
independent predictions, which can be made even in the absence of any information about 
the jet functions jW and the non-perturbative functions <pM, and without knowing 
A \onl M - For M = 7T°,p , we have A^ M = so Eq. fl^ gives 

For decays to D^K~, D^Kf', D^K°, and D^°Kf it was shown that AgJ = A£™ 



and so 



? k~ _ _ A(B° — > D*K~) 1 n xj- A(B°-Z>;*J 



-Rn — . ; ^„ „ S 77 x — 1 , Rn 



' A(B° -> D S K~) ' " -> A^f) 



The ratios in Eqs. (J5fij) and (JoTj) have calculable corrections of order a s (Q)/ir and power 
corrections 10 of order A/Q, which can be expected to be ~ 20%. 



10 Note that using the observed D and D* masses R(f = Nq/Nq = 1.04. This small difference corresponds 
to keeping an incomplete set of higher order corrections. 
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FIG. 5: The ratio of isospin amplitudes Rj = Ai/ 2 /{V2A 3 / 2 ) and strong phases 5 and <f> in 
B — > Dir and B — * D*ir. The central values following from the D and D* data in Table I are 
denoted by squares, and the shaded regions are the la ranges computed from the branching ratios. 
The overlap of the D and D* regions show that the two predictions embodied in Eq. H56|) work 
well. 



These amplitude relations imply the equality of the branching fractions. They also imply 
the equality of the non-perturbative strong phases between isospin amplitudes, namely the 
phases S D( * )M in the ratios Rf ( * >M as shown in Fig. |S| Thus for each of M — 7r°, p°, K°, K*® 

Br(B° - D*°M°) = Br(B° - D°M°) , 6 D '° M ° = 5 D ° M ° , (58) 

and for M = K~ , K*~ 

Br(B° -> D* S M) = Br(B° -> D S M) , 5 D ° M = 5 DsM . (59) 

The predictions in Eqs. (J56I58J) agree well with the data for D^n in Table HI which give 

|^| cx p = 0.94 ±0.21, 5 Dw = 30.3°^3 8 8 , 5 D * n = 30.1° ± 6.1° . (60) 

This agreement is represented graphically by the overlap of the la regions in Fig. with 
small squares indicating the central values. The dominant contribution to the phase 5 is 
generated by the (C — E) amplitudes which have complex phases from J® in Eq. (|55jh 
Since the phases in sf' 8 ^ are non-perturbative and can be large it is expected that they 
will dominate. Note that with this choice of triangle the power suppressed side in Fig. El is 
enlarged by a isospin pref actor of 3/ y/2 = 2.1. 

For B° decays to D^°p°, D^°K°, D^Kf, D^K' and D^Kf only upper bounds 
on the branching ratios exist, so our relation between D and D* triangles has not yet be 
tested. For each of these channels similar triangles to the one in Fig. El can be constructed 
once data becomes available. 

The results in Eqs. ()56j) and (|57|) can be contrasted with the absence of a definite pre- 
diction in the large N c limit as in Eq. (jSJ. Even when only the Fq term is included (naive 
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factorization), R n is given by a ratio of B — > 7r form factors, which for generic m^ c are 
not related by heavy quark symmetry. Thus, one does not expect a relation between the 
branching fractions or strong phases unless the 1/Q expansion is used. 

Next consider the second class of predictions, which follow from the perturbative expan- 
sion of the jet function in Eq. (|55|). We now assume that a g (fi ) is perturbative, and focus 
on M = it, p since the kaons are contaminated by contributions from A^^ M . The tree level 
result for J is given in Eq. (J42)) . and was used to define the nonperturbative parameters 
s (o,8) through convolutions with the soft distribution functions S^' 8 \kf) as in Eq. (|44|) . It 
is convenient to introduce an effective moment parameter, 

i /~»(8) 

s „ = _ s (o) _j_ t r^ s W = \ s _| e -tf (61) 

SeS S + 2NcCfc{0 ) S l S eff|e • IDlj 

L 

In terms of the effective moment the result in Eq. (J55)l at lowest order in a s (Q) and a s (po) 
becomes 

A(B° - DWM°) = N M cf 16TO n s(/io) Seff ( M o) {x- X )m , (62) 

where N^ 1 is defined in Eq. Since s c s is independent of M = ir, p the same phase <fi is 

predicted for these two light mesons. 

At leading order in 1/Q the type-I amplitude Ao~ = A(B~ — > D n~) factors as in Eq. (J7J) 
giving the product of a form factor and decay constant, both of which are real (with the 
usual phase conventions for the states, and neglecting tiny a s (rrib) strong phases (~ 2°) 
generated by the coefficients Cf^ R at one-loop Q). Therefore the amplitude A _ is real at 
leading order in l/Q, up to calculable corrections of order a s (Q). Choosing the orientation 
of the triangle so that A _ lies on the real axis, the phase <ft can be directly extracted as one 
of the angles in the isospin triangle 

V2A 00 + A + _ = A _ . (63) 

This is shown in Fig. |U1 where we divide by A _ to normalize the base. The data on B° — > 
D°p° is not yet sensitive enough to test the prediction that <fi is the same for 7r° and p°. 

Using Eqs. ((7j) and (jU2*|) it is possible to make a prediction for the ratio R c in Eq. (|12|) at 
NLO in the power expansion. Since R c = A + _/A _ contains only charged light mesons it 
is easier to measure than neutral pion channels. Data is available for all four of the D^*\ 
and D^p channels. Using the triangle relation in Eq. (JBlIj) one finds for the ratio of any two 
such modes [M = ir,p] 

A _ 9E M (m B + m DM ) Z(w , p ) V /M 1 ; 

It is easy to see that the ratio of amplitudes on the right-hand side is common to final 
states containing aDor D*, and has only a mild dependence on the light meson, introduced 
through the inverse moment (x -1 )^/. In particular we note that there is no dependence 
on the decay constant fu on the RHS of Eq. (JEU), since it cancels in the ratio A 00 /A -. 
This implies that the ratios R c are comparable for all four channels D^tc and D^p, up 
to corrections introduced by (x^ 1 )^ ^ (x~ 1 ) p . These corrections can be smaller than the 
correction one might expect from the ratio of decay constants f p /f n — 1.6 (which appear in 
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the naive factorization). The experimental values of these ratios can be extracted from 
Table I and are in good agreement with a quasi-universal prediction 



\R{ Dn) \ 

\Ri D ^\ 
\R™\ 



\R 



( d IPl)\ 



\A(B° - 


/i 


\A(B~ - 


-> D°7r- n 


\A(B° - 


-> D* + TT-)\ 


\A{B~ - 


^D*°n-)\ 


\A(B° - 


+ D*p-)\ 


\A(B~ - 


-Z>V)| 


\A[B° - 


-D*+p~)\ 


\A{B~ - 


-> D*°p~)\ 



0.77 ±0.05 



0.81 ±0.05, 



0.80 ±0.09, 



0.86 ±0.10. 



(65) 



This lends support to our prediction for the universality of the strong phase <fi in B — ► D^n 
and B — > D^p decays from the s e s in Eq. (|64j) . The central values of R c ~ 0.8 are well 
described by an s e fj of the expected size (~ Aq C d) as discussed in the fit to the isospin 
triangle below. Further data on these channels may expose other interesting questions, such 

as whether Ri° M ^ is closer to R^ M>> than r[ d( ^ is to R[ D< )p \ 

An alternative use of Eq. (|64|) and the R c amplitude ratios is to give us a method for 
extracting the ratio of p and 7r moments. Using the Dn and Dp measurements which have 
smaller errors than for D*, we find 



\R { c Dp) \ 



1 



\R ( c D7r) \ 



0.87 ±0.42. 



(66) 



where only the experimental uncertainty is shown. The extraction in Eq. ()66)1 is smaller, 
but still in agreement with the ratio extracted from light-cone QCD sum rules. The best fit 
from the 7*7 — > ir° data performed in Ref. gives (x -1 ),,- = 3.2 ± 0.4 in agreement with 



sum rule estimates of the moment. The QCD sum-rule result (x l ) p = 3.48 ±0.27 |40(, then 
implies 



1.10 ±0.16. 



(67) 



SR 



The result that this ratio is close to unity is consistent with the universality of the data in 
Eq. (jnHJ). This data can be contrasted with cases where the single light meson is replaced 
by a multibody state such as 



Br(B c 



D 



7T 7T 7T 7T 



0> 



Br(B- 



1.02 ±0.27. 



7T 7T 7T L 



(68) 



For the four pion final state our proof of the factorization theorem does not work, since 
for many events one or more of the pions will be slow. We therefore would expect less 
universality in branching ratios involving more than one light meson. (For these decays a 
different type of factorization involving large N c works well for the q 2 spectrum |8|.) 

The result in Eq. (Jo^j) also leads to predictions for the ratios of color-suppressed decay 
amplitudes to final states containing different light mesons M° = vr°,p°. We find 



dp/ 71 " — 

it n = 



\A(B°^D°p°)\ _ f p (x- 1 ), 



1.40 ±0.77. 



(69) 
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where we used f n ± = 130.7 ± 0.4 MeV, and f p ± = 210 ± 10 MeV, and inserted the result in 
Eq. (jUBj) for the moments. This can be compared with the experimental result (R p / n ) cxp = 
1.02 ± 0.21. The large uncertainty in the ratio of moments in Eq. ()66|) dominates the error 
in Eq. With the QCD sum rule result in Eq. (jUZD we find R p/w = 1.64 ± 0.35, a result 
whose central value is farther from the experimental data, but still consistent with it. 

In contrast to the first class of predictions, the predictions for the ratios in Eqs. (|64|) . (|66|). 
and (|69|) and the prediction for the universality of can receive corrections from neglected 
[a s (/io) 2 /ii] terms in J. The dominant theoretical corrections to this extraction are expected 
to come again from these perturbative corrections to J or from power corrections, which 
we estimate may be at the ~ 30% level. A future study of the perturbative corrections is 
possible within the framework of our factorization theorem and SCET. 

The result in Eq. ()62|) and the data on B — > Dtt and B —>■ D*tt decays can be used to 
extract values of the moment parameters |s e ff | and strong phase 0. We present in Fig. El 
the constraints on the parameter s e g in the complex plane, obtained from Dtt (light shaded 
region) and D*tt data (darker shaded area). We used in this determination fi = E w = 
2.31 GeV, and leading order running which gives a s (/i ) = 0.25, Ci(fi = /i ) = 1.15, and 
Ci{\i = A*o) — —0.32. The good agreement between the Dtt and D*tt la regions marks a 
quantitative success of our factorization relation in Eq. (|55|). Averaging over the Dtt and 
D*tt results, we find the following values for the soft parameters at // = //q 



|s eff | = (428 ± 48 ± 100 MeV 



0.26 \ / 3.2 



C L (fj, )a a (fj, )J \(x 1 } 7T/ 
= 44.0° ±6.5°. (70) 

In this determination the inverse moment of the pion wave function was taken from the best 
fit to the 7*7 — > 7T° data [39], (a:™ 1 )^ = 3.2±0.4. For \s e g | the first error is experimental, while 
the second is our estimate of the theoretical uncertainty in the extraction from varying /j,q 
from En/ 2 to 2E n . At the order we are working the extraction of the phase is independent 
of the scale, since the prefactor a s (fi )(x^ 1 ) 7T drops out. The result in Eq. (|70j) agrees well 
with the dimensional analysis estimates s e s ~ s^ ' 8 -* ~ Aqcd- Since is non-perturbative its 
value is unconstrained, and a large value of this phase is allowed. 

The recent B° — > D°p° data from Belle allows us to extract |s e ff| and in a manner 
independent of the above determination. Keeping only experimental errors we find 



|seff| = (259 ±124 MeV 



0.26 \ / 3.5 



C L (fi )a s (fi ) J \{x !) f 
= 17° ±70°. (71) 

The results agree with Eq. ((70)) within la, but currently have errors that are too large to sig- 
nificantly test the factorization prediction of equality on the 20-30% level of the parameters 
extracted from Dp and Dtt. 

The B° — > Ds K~ channels proceed exclusively through the W^-exchange graph and 
have been the object of recent theoretical work 44]. For the result analogous to Eq. (|62)1 
we would have [M = K, K*} 

A(B° - D^M) = V2N M Cf l ±^AtA s e M {x -i )m + j^m (72) 
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FIG. 6: Fit to the soft parameter s e g defined in the text, represented in the complex plane 
with the convention that Aq_ is real. The regions are derived by scanning the la errors on the 
branching fractions (which may slightly overestimate the uncertainty). The light grey area gives 
the constraint from B — > Dtt and the dark grey area gives the constraint from B — > D*tt. 

Both the B° — > K~ modes and the Cabibbo-suppressed decays B — > _D(*)j^(*) receive 
this additional contribution from A^^ M . This makes the factorization theorem less pre- 
dictive, and so we do not attempt an analysis of ratios R^* )K{ *\ R^ M , or the universal 
phases 4>e and (pc that are analogous to the cf> in Eq. (jfTTj) . 

On the experimental side both the Babar and Belle Collaborations 0] recently observed 
the B° — > D^K~ decay, and set an upper limit on the branching ratio of B° — > D* + K~ 

B(B° -> D+K-) = [3.2 ± 1.0 (stat) ± 1.0 (sys) ] x 10~ 5 (Babar) (73) 

= [4.6t^ (stat) ± 1.3 (sys) ] x 10" 5 (Belle) 
B(B° -> D* s + K-) < 2.5 x 10" 5 (90%CL) (Babar) . 

The branching fraction for B° — > D+K~ is an order of magnitude smaller than that for 
B° — > D 7T°. This indicates that the ^-exchange amplitude E DsK is suppressed relative to 
(C - E) Dn and (V ud /V2V us ) C D ° k ° . In SCET the SU(3) breaking between <j> v {x) and <f> K (x) 
is generated by masses in the collinear quark Lagrangian ji^. This causes an asymmetry 
in the light-cone kaon wavefunction. This SU(3) violation can be expected to be at most a 
canonical ~ 20-30% effect, which would not account for the observed suppression. 

However, there is one important source of potentially larger SU(3) breaking from an 
enhancement in moments of the light-cone kaon wavefunction which appear in the short 
distance amplitude. This may account for the observed suppression. Basically strange 
quark mass effects imply a larger SU(3) violation for inverse moments than expected for 
07r versus <pK alone, and implies that {x^)k < { x ^ 1 )k- Using the result from QCD sum 
rules the ratio of moments is (x^ 1 ) k / X^J 1 ) k ~ 1-4. Furthermore, we anticipate a 
similar large effect from the moments that appear in the soft matrix elements which again 
differ by factors of (fcj~) -1 versus (fc+) _1 , and appear in a way that suppresses D S K~ . The 
combination of these two suppression factors might accommodate the observed factor of 



rn= Dtt 
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three suppression in the D S K~ amplitudes. 11 The long distance amplitude also involves 
two inverse momentum fractions through j^ 0,8 ^ i n Eq. (j2HJ), although admittedly much less 
is known about the non-perturbative functions &m and Q^'r '■ Thus, we find that the 
suppression of E DaK may not imply much about the relative size of C Dn and E Dn . Finally, 
we note that the suppression mechanism for ss creation that we have identified is particular 
to problems involving large energies where light-cone wavefunctions arise. 

Further information on the relative size of the short and long distance contributions to 
the kaon factorization theorem is clearly desirable. In section ITVl it was noted that in type-II 
decays transverse K*'s are produced only by the long distance contribution at this order 
in Aqcd/Q- Therefore, measuring the polarization of the K* in both the B° — > D*K*~ 
and B° — > D*°K*° decays can give us a direct handle on whether there might be additional 
dynamical suppression of either the long or short distance contributions, or whether they 
are similar in size as one might expect apriori from the power counting. 



VII. SUMMARY AND CONCLUSIONS 



We presented in this paper the first model-independent analysis of color-suppressed B° — ► 
_D°(*)M° decays, in the limit of an energetic light meson M°. The soft-collinear effective 
theory (SCET) was used to prove a factorization theorem for these decay amplitudes at 
leading order in Aqcd/<5, where Q = {mj, m c , Em}- Compared with decays into a charged 
pion these decays are suppressed by a factor Aqqb/Q- Therefore, in the effective theory 
they are produced exclusively by subleading operators. 12 

We have identified the complete set of subleading operators which contribute to B° — > 
£>o(*) M o decays with M = vr, p, K, K*, as well as for the decays B° -> D { * ] . After hard- 
soft-collinear factorization, their matrix elements are given by i) a short distance contribution 
involving a jet function convoluted with nonperturbative soft distribution functions, and 
the non-perturbative light-cone meson wave function, and ii) a long distance contribution 
involving another jet function and additional x± dependent nonperturbative functions for 
the soft B, D and collinear M. The long distance contributions were shown to vanish for 
M = ii, p at lowest order in a s (Q)/ir. 

The factorization formula is given in Eqs. and (j37|) . It may seem surprising that the 
type-II decays factor into a pion light-cone wave function and a B — > £)W soft distribution 
function rather than being like the naive factorization in Eq. (j2J). Our results indicate 
that factorization for type-II decays is similar to factorization for type-I decays (albeit with 
new non-perturbative soft functions and additional long distance contributions for kaons). 
To derive Eq. (|35|). QCD was first matched onto SCETj at the scale p 2 = Q 2 . In SCETj it 



In general this argument gives a dynamic explanation for the suppression of ss-popping at large energies 
which could be tested elsewhere. The production of an ss pair which end up in different strange hadrons 
is likely to be accompanied by a suppression from inverse momentum fractions that arise from the gluon 
propagator that produced these quarks. This enhances the SU(3) violation in a well defined direction so 
that less ss pairs are produced. A factor of 3 suppression of ss popping is implemented in JETSET [4^ . 
In type-I decays, other subleading operators can compete with the time ordered products we have identified 
at the same order in A/Q. This makes a complete analysis of power corrections to type-I decays more 
complicated than our analysis of type-II decays. 
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is still possible for gluons to redistribute the quarks. This intermediate theory provides a 
mechanism for connecting the soft spectator quark in the B to a quark in the pion, and for 
connecting the energetic quark produced by the four-quark operator with the soft spectator 
in the D (see Fig.|2J). This process is achieved by the power suppressed time ordered products 
given in Eq. (|21|). SCETj is then matched onto SCET n at a scale /1q = E M A. In SCET n the 
collinear quarks and gluons are non-perturbative and bind together to make the light meson 
M. This second stage of matching introduces a new coefficient function (jet functions) as 
in Eq. (|24j) . The jet function J contains the information about the SCET I graphs that 
move the spectator quarks into the pion. The physics at various scales is neatly encoded 
in Eq. ([35)1 . The Wilson coefficient C(z) from matching QCD onto SCET I depends on 
physics at the scale Q 2 , the jet functions J, J from matching SCET I onto SCET II depends 
on QA physics which is where quark redistribution occurs, and finally the soft distribution 
functions S, $ and the pion light cone wavefunction 4>m,^m depend on non-perturbative 
physics at A 2 which is where the binding of hadrons occur. 

The soft functions 5* are complex, and encode information about strong rescattering 
phases. This information is introduced through Wilson lines along the light meson direction 
of motion, which exchange soft gluons with the final state meson D^*\ They provide a 
new mechanism which generates non-perturbative strong phases. In the literature other 
mechanisms which generate perturbative strong phases have been proposed. In particular 
in Ref. 0, H3| a method for identifying perturbative strong phases with an expansion in 
a s (Q 2 ) was developed. In Ref. |38|, |46| it was pointed out that strong phases can also 
be generated perturbatively at the intermediate scale a s (E M A). In the language of our 
factorization theorem in Eq. (J35)) these phases roughly correspond to imaginary parts in 

the hard coefficients C^'^ and jet functions J respectively. These phases exist, but for 
the B — > Dit channels they only show up at next-to-leading order in the a s {m b ) or a s {uo) 
expansion. (In type-I B — > D^n decays the hard strong phase is very small, ~ 2° |7[). 
In contrast, our new source of strong phases is entirely non-perturbative in origin and can 
produce unconstrained phases. For the case of B — > D^M these phases show up in the 
power suppressed class-II amplitudes. 

The factorization theorem proven in this paper leads to predictions which were tested 
against existing experimental data on color-suppressed decays. We derived two model inde- 
pendent relations, which related 

• the B° — > D°M° and B° — > D*°M° decay branching fractions and 

• the B — > DM and B — > D*M strong phases. 

Here M = tt, p, K, Kp and these relations are true to all orders in the strong coupling at 

the collinear scale. The same predictions are also obtained for B° — > Dg K~ and B° — > 
D s K,, . The good numerical agreement observed between the strong phases and branching 
fractions in the Dir and D*ir channels gives strong backing to our results. This prediction 
can be tested further since the equality of the strong phases for the p, K, and Kt channels 
have not yet been tested experimentally. 

Additional predictions followed from the factorization theorem by using a perturbative 
expansion for the jet function, including [M = n, p] 

• the ratios 1^1 = \A(B° -> D^ + M-)/A(B- -> D^°M~)\ to subleading order 

• the ratios \R p /n \ = \A(B° -> D^°p°)/A(B° -> D^°tt°)\ to subleading order 
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• universal parameters {|s e fj|,0} which appear for both D^n and D^p, and 

• a mechanism for enhanced SU(3) violation in ss production for the short distance 
amplitude which might explain the suppression of the B° — > DgK~ rates relative to 

For \R C \ taking different values of M with the same isospin the power corrections only 
differ by the moments (x~ 1 )m, giving an explanation for the observed quasi-universality of 
these ratios. The isospin triangles for these M's are predicted to involve a universal angle 
0. The ratio of neutral modes I-Rq^I are determined by inverse moments of the light-cone 
wavefunctions and decay constants. Finally extractions of the non-perturbative soft moment 
parameter s e s agrees with the ~ Aq C d size estimated by dimensional analysis. 

In the case of B° D [ * ] K^- an additional suppression mechanism was identified, which 
arises from enhanced SU(3) violation due to the asymmetry of non-perturbative distributions 
involving strange versus down quarks. The inverse moments that appear in the factorization 
theorem enhances this difference, and can lead to a dynamic suppression of ss-popping. 
Further information on the size of the short and long distance amplitudes would help in 
clarifying this observation. 

A more detailed experimental study of the channels in Tables I and II is crucial to 
further test the accuracy of the factorization theorem and improve our understanding of 
the structure of power corrections. Work on extending these results to decays to isosinglet 
mesons is in progress. It should be evident that B s decays could also be considered although 
we have not done so here. 
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APPENDIX A: LONG DISTANCE CONTRIBUTIONS FOR vr AND p 

The factorization theorem derived in Sec. lIIII for the color-suppressed B° — > D°M° ampli- 
tude contains both short- and long-distance contributions. In this Appendix we show that, 
working at lowest order in the Wilson coefficients at the hard scale Q, the long-distance 
amplitude vanishes for the case of an isotriplet light mesons M = tc, p. 

We start by recalling the factorized form of the long-distance amplitude, which is given 
by SCET n time ordered products T^'j^ 

Aong M = J o dzJdk + du;Jd 2 x ± [c£\z) JW(w*+) ^2(^^,4/) 

±Cf (z) P\ook + ) $%(k+x x ,e^)*®(z,L;,x x ,e* M )\ . (Al) 

The functions and are SCETn matrix elements of collinear and soft fields, re- 

spectively, and their precise definitions are given in Eqs. flUEJ). The jet functions J^(uk + ) 
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appear in the definition of the subleading soft-collinear Lagrangian and their lowest 

order expressions are given in Eq. ([26)1 . 



In the following we derive a few general properties of the functions 



M 



and 7 ( °' 8) fol- 

lowing from isospin, charge conjugation, parity and time-reversal. The collinear function 
*&m(z, u, s* m ) is defined as the matrix element 



M°(e) 



[&w) t 4p l {w^^) T2 ] (o ± ) \&w)^p L (w^^\] ( X± ) 



(A2) 



We will prove that ^ M=7T p is even under uo — > —uo and z — > 1 — z. As motivation consider 
the first bilinear in Eq. (jA2|) . which creates a du collinear quark pair. The second bilinear 
in Eq. (|A2|) must act at some point along the collinear quark lines: it either takes a d — > u 
(for uj > 0) or takes a w — > (for u; < 0). Examination of lowest order graphs contributing 
to \I/m shows that these two types of contributions always appear in pairs, such that the 
projection of \1/a/ onto an isotriplet state is even under uo — > —uo. This suggests the existence 
of a symmetry argument, valid to all orders in perturbation theory. 

We will prove that *£' 8) is even, as a consequence of G-parity. This is defined as usual 
by G = C exp(— inh) where C is charge conjugation and I2 is the isospin generator, and is a 
symmetry of the collinear Lagrangian in the limit m u ^ <C Aqcd- Its action on the collinear 
operators in Eq. (jA2j) can be worked out from that of its components C and J 2 (cf. Ref. [28j]) 
and is given by 



(d)\ 



(A3) 



Taking into account the G-parity of the states, Eq. (|A2|) is equal to 



±(M°(e 



(^WUJPr(W^^)J (Oj.) \&W)^P R {WH^)-J\ (x±) 0) , (A4) 



where the ± refer to the p° and 7r° respectively. Next we apply parity in the matrix element 
followed by switching our basis vectors n «-> n. Acting on Eq. (jA4j) this gives 



&W)_ T ^P L (W^^)„ T1 (0 ± ) &W)^$P L {W^. 



-x±. 



(A5) 



where the overall sign is now the same for M = p,n. 
the only allowed perpendicular dot products are (—x^ 2 



Now since \E^ 8 ^ is a scalar function 



x]_ and — x±-e* P = x^-e* . Finally 



we note that the change in n 2 from Eqs. (jA2j) to (jA5|) is equivalent z — > 1 — z. Thus the 
invariance of SCET n under G-parity and regular parity has allowed us to prove that 



TT.p V 



*® p (l-z, -u,x ± ,e*). 



(A6) 



-7(0,8) ( 



Next we prove that J v "'" (co»A; + ) is odd under uj —>■ —uj. By reparameterization invariance 
type-Ill |47| only the product uok + will appear. Consider applying time reversal plus the 
interchange (n <->• n) to the SCETn Lagrangian. Since this Lagrangian does not have coef- 
ficients that encode decays to highly virtual offshell states it should be invariant under this 
transformation. Acting on Eq. (|25|) this implies that J^ ' 8 ^ must be real, 



[j m (uk + )y 



-7(0,8), 

J [uok 



(uuk + ). 



(A7) 
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At tree level this implies that we should drop the ie in the collinear gluon propagator in 
matching onto this operator. This was done in arriving at the odd functions j( 0,8 ) oc l/(cok + ) 
in Eq. (|26|). The imaginary part would give a 5(uk + ) and corresponds to cases where the 
SCETx T-product is reproduced by a purely collinear SCET n T-product {k + = 0), or a 
purely soft SCET n T-product {uj = 0). Thus dropping the ie also saves us from double 
counting. 

Now consider what functions can be generated by computing loop corrections to t / ' 8 ' 1 . By 
dimensional analysis t / ' 8 ' ) must be proportional to l/(uk + ) times a dimensionless function 
of uk + /fi 2 . Since at any order in perturbation theory the matching calculation will involve 
only massless quarks we can only generate logarithms. Therefore, we must study functions 
of the form 



To demand that only the real part of these functions match onto J ' we average them with 
their conjugates. It is straightforward to check that only terms odd in uj — > — u survive. 

Thus, all the terms that can correct the form of J^°' 8 at higher orders in a s are odd under 
uj — > — to. 

Now in Eq. (|A1|) the integration over uj is from — oo to oo, while z varies from to 1. 
Consider the change of variable uj — > — uj and z — > 1 — z. If C^ R (z) = C^ R (1 — z) then 
under this interchange one of the functions in the integrand is odd (J) and the other two 
are even (C^ R and so the integral would vanish. 

Now if Cf R { z) are kept only to leading order then they are independent of z and thus 

unchanged under z —> 1 — z. So at this order in the a s (Q)/ir expansion of C^' R (z) we 

find A^ g M = 0. This completes the proof of the assertion about the vanishing of the long 
distance contributions for M = 7T, p. 

APPENDIX B: HELICITY SYMMETRY AND JET FUNCTIONS 

In this appendix we discuss the general structure of the jet functions J(°> 8 \z, x, k^) in 
Eq. (|35[). which are generated by matching SCETj and SCETn & t any order in a s (fi ). In 
Fig. |2k,b this means adding additional collinear gluons which generate loops by attaching 
to the collinear lines already present (as well as vacuum polarization type collinear quark, 
gluon, and ghost loops). Additional collinear loops should also be added to Figs. Et,d,e, 
and the difference at lowest order in A gives J. Throughout this appendix we continue to 
drop isosinglet combinations of £ n • • • £ n . These will also have additional contributions from 
topologies where the outgoing collinear quarks are replaced by outgoing gluons (through 
operators). 

The leading order collinear Lagrangian has a U(l) helicity spin symmetry for the quarks, 
see reference 2 in j2^|. It is defined by a generator h n that has the quark spin projection 
along the n direction, which is different from usual definition of helicity as the projection of 
the spin along its momentum. Unlike QCD, the collinear fields in SCET only allow quarks 
and antiquarks that move in the n direction. For h n we have 



&» = Te±<V, h 2 n = l, [h n ,i] = [h n ,ifi\ = 0, {h n ,Y±} = 0- (Bl) 




(A8) 



1 
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After making a field redefinition |26( to decouple ultrasoft gluons the leading order collinear 
quark Lagrangian is 



4? = W \ in ■ Dc+^t-^rr^ \ ? (B2) 



where iD% contains only collinear A*^ q gluons. is invariant under the transformation 

£ n — > exj>(i8h n )£ n , £ n — > £ n exp(—i9h n ). This means that any number of leading order 
collinear quark interactions preserve the quark helicity h n . The collinear gluon interactions 
take w n (t) -> Mn(T), Wnd)-^ Ml)> v n {\) -> u n (t)> u n(j) -> ^n(l), and can also produce 
or annihilate the quark-antiquark combinations w n (T) u n(l) or w n (|) f n (t) (the arrows refer 
to the helicity of the antiparticles themselves rather than their spinors). For this reason we 
refer to clS 8b Ah n = operator. 

The leading order SCETi operators in Eq. (j!8|) are also unchanged by the h n - 
transformation and therefore does not change collinear quark helicity. In contrast the op- 
erators CP do generate or annihilate a collinear quark giving Ah n = ±1/2. However, at 
tree level we showed in Section 11111 that the two graphs in Figs. |2K,b match onto an overall 
Ah n = operator in SCET n as given in Eqs. (|25|) . Since at higher orders the C^) will not 
cause a change in the helicity they also match onto these same operators, so the structure 
will not occur. At tree level only the structure ?/iPl ® %>Pl appeared in Eq. (|28|). To 
rule out the appearance of Pr beyond tree level we note that the weak operator projects 
onto left handed collinear fermions, and for the jet function the conservation of helicity in 
Cgj: implies a conservation of chirality. This leaves us with the desired result. 

It is perhaps illustrative to see this more explicitly by looking at the spin structure of 
the loop graphs. We begin by noting that the spin and color structure in h$ ■ ■ ■ hv is 
unaffected by this second stage of matching. Adding additional collinear attachments only 
can affect the spin and color structure generated in putting the collinear quark fields and 
light ultrasoft quark fields together. 

Consider how additional gluon attachments effect the spin structures that appear in 
Figs. |2^,b. The leading order collinear quark Lagrangian is Cft? in Eq. (JB2J). Each attach- 
ment of a collinear gluon to a collinear quark lines in the figures generates a ^/2 from the 
vertex and a ifij2 from the quark propagator. These combine to a projector which can be 
eliminated by commuting them to the right or left to act on the collinear quark spinors, via 
(^)/4 £ n = £ n . Therefore, at most we have additional pairs of 7^'s that appear between the 
light quark spinors. The aim is to show that just like the tree level calculation in Eq. 
the resulting operators have spin structure (difiPiu) (£ n ^-Pz£n)- 

For the contraction of Tj' which gives the C topology the spin structure is 

u^iTiT ■ ■ - t^tT (I^)tI 1 • ■ ■ iT (f iiy s u) ] (71)7^ • ■ ■ 7Nf 

= [u { n d hT ■ ■ ■ ll 2k 'llP L u^] [u^hhl 1 ■ ■ ■ 7Ni d) ] • (B3) 

In the first line (f-Pz) comes from Q^r , the 7° and 7^ are terms generated by the 
insertions and the ifi/2 is from the extra collinear quark propagator. In the second line the 

Pl projector was moved next to Us without a change of sign (for anticommuting 75), and 

the remaining t/i and ijt were then moved next to the u$ and canceled. The remaining free 
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_L indices in the second line are contracted with each other in some manner. Fierzing the 
set of 7 matrices in Eq. (jB3|) by inserting 1 ® 1 next to the collinear spinors gives 

H d) ri v<£>] H'hhi 1 ■ ■ ■ 7^ ri iT ■ ■ ■ i^'ilP^} , (B4) 

where 

ri®r; = | ®i- ®^ 75 _ M ^ 

- |(l-7fe)®lt-^®^. (B5) 

In the second line of Eq. (jB5|) we have used the fact that the 75 in the bracket with soft 
quark spinors can be eliminated by moving it next to the P^. To eliminate the ^7^ Dirac 
structure we note that between the soft spinors in Eq. (jB4|) there are an odd number of 
7_i_'s to the left and right of ^7^, and so at least one set of indices are contracted between 
the sets {/3, Ai, . . . , Xzj} and . . . , //2fc', The identity {7x ; 7i} — can be used 
to move these matrices so that they sandwich 7^, and this gives the product 7^_7^7^~ = 0. 
After these manipulations only the spin structure (di^Piu) {C.n'^PL^n) remains. A similar 
argument can be applied to the E-topology with the same result. 

In several places in the above argument we made use of Dirac algebra that is particular 
to 4-dimensions (anticommuting 75 and setting 7^7,^7^" = 0). If the 7_i_'s are taken in in full 
dimensional regulation then it is not apriori clear if the manipulations survive regulation. 
However, the original helicity symmetry argument shows that as long as the theory can be 
regulated in a way that preserves this symmetry this will indeed be the case. 



APPENDIX C: PROPERTIES OF SOFT DISTRIBUTION FUNCTIONS 

In this appendix we derive some useful properties of the soft functions S^ 0,8 ) . In particular 
we show that these functions are complex. The imaginary parts have a direct interpreta- 
tion as non-perturbative contributions to final state rescattering between the and final 
energetic meson as discussed in section 11111 

To be definite we consider the function \ and suppress the index L. The manipulations 
for the remaining soft functions and S^ R are identical. The definition in Eq. (|34|) is 

(D°(^OI(^ ) 5)^P i (5t/ i f)(JS') fci+ ^(5'M A2+ |S ^)^ = S<°>, (CI) 

where the Wilson lines are defined as 

W = [ ^ exp ( - £ n-A n , q (x) )] , S = [ ]T exp ( -g J- n-A s ^] . (C2) 



perms perms 



In general is a dimensionless function of v-v', n-v, n-v', n-ki, n-k 2 , Aq CD , and \i. 
Since (S^q) k + = 5(k£ — n-V)(S*q) the LHS is invariant under a type-Ill reparameterization 
transformation [47] (n — > e a n, n — > e~ a n). Therefore the RHS can only be a function of w, 
t = n-v/n-v', z = n-ki/n-k2, K//i = [n-k\ n-fc 2 /(n-w n-v'/j 2 )} 1 ^ 2 , and Aqcd/a^- 
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Rather than study the matrix element in Eq. (|C1|) directly it is useful to instead consider 
(^(^)|(^5)^P L (St^)(gS) fc ^P L r a (Stg) fc2+ |^(t;)) 

= (G3) 

where h v are doublet fields under heavy quark flavor symmetry, and q and \H i= i^{v )) are 
isospin doublets of (u, d). The last three variables in Eq. (jC3|) will not play a crucial role so 
we will suppress this dependence. Taking the complex conjugate of Eq. (|C3J) gives 

= (C4) 

The dependence on u> and K is unchanged since they are even under the interchange v <-> t>', 
n-fci <-> n-/c2- Next, decompose the functions in terms of even and odd functions under 

t -> 1/t, 2 -> 1/^: 

5 (o) = 5 g) + 5 (o) ^ (C5) 
where S { ^o = [S (0) (M) ± S<®(l/t, l/z)]/2. Now Eq. Jd implies that 

[fif Ml* = fifM, = -4 0) (*^) (C6) 

so S 1 ^ is real and S'q is imaginary. An identical argument for implies that it too is a 
complex function. 

For the above analysis it is important to note that n ■ v ' = mB/mn is not 1 in the heavy 
quark limit where we have new spin and flavor symmetries. These symmetries arise from 
taking m B ^> Aqcd and mn ^> A QCD , not from having mg = m D . 
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